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DRIFT- FREE VIDEO ENCODING AND DECODING METHOD, AND CORRESPONDING DEVICES 



FIELD OF THE INVENTION 

The present invention relates to an encoding method for the compression of an 
original video sequence divided into successive groups of frames (GOFs) and to a 
corresponding decoding method. It also relates to corresponding encoding and decoding 
devices. 

BACKGROUND OF THE INVENTION 

The growth of the Internet and advances in multimedia technologies have 
enabled new applications and services. Many of them not only require coding efficiency but 
also enhanced functionality and flexibility in order to adapt to varying network conditions 
and terminal capabilities. Scalability answers these needs. Current video compression 
standards often use so-called hybrid solutions, based on a predictive scheme where each 
frame is temporally predicted from a reference frame (the prediction options being : zero 
value prediction, for the intra frames or I frames, forward prediction, for the P frames, or bi- 
directional prediction, for the B frames) and the obtained prediction error is spatially 
transformed to get advantage of spatial redundancies. From MPEG-2 to MPEG-4, standard- 
based scalable solutions have then been proposed. They rely on the generation of a base 
layer, containing the lowest spatial, temporal and/or SNR resolution version of the original 
video sequence, and one or several enhancement layers allowing (if transmitted and decoded) 
a spatially, temporally and/or SNR refined reconstruction. A short-coming of these layer- 
based scalability schemes comes however from their lack of coding efficiency. 

A different approach has been proposed with techniques such as three- 
dimensional (3D) subband coding, which are able to generate embedded bitstreams. Thanks 
to their multi-resolution analysis structure, scalability is inherent to these schemes and does 
not weaken their intrinsic coding efficiency. In a 3D subband codec such as described for 
example in "A fully scalable 3D subband video codec", "Proceedings of the International 
Conference on Image Processing (ICIP2001), vol.2, 2001, pp.1017-1020, the embedded 
bitstream is fully scalable and can be decoded at any spatial and temporal resolutions, and 
with any desired SNR quality, simply by truncation at known locations. In such a scheme, 
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successive groups of frames (GOFs) are processed as a 3D structures and spatio-temporally 
filtered in order to compact the energy in the low frequencies, a motion compensation being 
also provided in order to improve the overall coding efficiency. The 3D subband structure is 
depicted in Fig.l : the illustrated 3D wavelet decomposition with motion compensation is 
5 applied to a group of frames (GOF), and this current GOF is first motion-compensated (MC), 
in order to process sequences With large motion, and then temporally filtered (TF) using Haar 
wavelets (the dotted arrows correspond to a high-pass temporal filtering, while the other ones 
correspond to a low-pass temporal filtering). After the motion compensation operation and 
the temporal filtering operation, each temporal subband is spatially decomposed into a 

10 spatio-temporal subband, which finally leads to a 3D wavelet representation of the original 
GOF, three stages of decomposition being shown in the example of Fig.l (L and H = first 
stage ; LL and LH = second stage ; LLL and LLH = third stage) .The well known SPIHT 
algorithm, extended from 2D to 3D, is chosen in order to efficiently encode the final 
coefficient bit-planes with respect to the spatio-temporal decomposition structure. 

15 As it is implemented now, a 3D subband codec applies the 

motion-compensated (MC) spatio-temporal analysis at the full original resolution at the 
encoder side, Spatial scalability is achieved by getting rid of the highest spatial subbands of 
the decomposition. However, when motion compensation is used in the 3D analysis scheme, 
this method does not allow a perfect reconstruction of the video sequence at lower resolution, 

20 even at very high bit-rates : this phenomena, referred to as drift in the following description, 
lowers the visual quality of the scalable solution compared to a direct encoding at the 
targeted final display size. As explained in the document 'TVIultiscale video compression 
using wavelet transform and motion compensation", P.Y.Cheng and al., Proceedings of the 
International Conference on Image Processing (ICIP95), Vol.l, 1995, pp.606-609, this drift 

25 comes from the order of wavelet transform and motion compensation that is not 

interchangeable. Indeed, when a frame (A) is synthesized at a lower resolution (a), the 
following operation is applied : 



30 where DWT L denotes the resolution downsample using the same wavelet filters as in the 3D 
analysis. In a perfect scalable solution, one wants to have: 



The remaining part of the expression (1) therefore corresponds to the drift. It can be noticed 
that, if no MC is applied, the drift is removed. The same phenomena happens (except at the 



a = DWT L (L)+ MC[DWT l (H)] 

= DWT L (A)+ [MC[DWT L (H)] - DWT L (MC[H])] 



(1) 



a = DWTl (A) 



(2) 
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image borders) if a unique motion vector is applied to the frame. Yet, it is known that MC is 
unavoidable to achieve a good coding efficiency, and the likelihood of a unique global 
motion is small enough to eliminate this particular case in the following paragraphs. 

Some authors, such as J.W.Woods and al in the document "A resolution and 
5 frame-rate scalable subband/wavelet video coder", IEEE Transactions on Circuits and 
Systems for Video Technology, vol.l, n°9, September 2001, pp.1035-1044, get rid of this 
drift to achieve good spatial scalability by different means. However, in said document, the 
described scheme, in addition to being quite complex, implies the sending of an extra 
information (the drift correction necessary to correctly synthesize the upper resolution) in the 
10 bitstream, thus wasting some bits (the solution described in the document "Multiscale video 
compression. . avoids this bottleneck but works on a predictive scheme and is not 
transposable to the 3D subband codec). 



SUMMARY OF THE INVENTION 
15 It is therefore an object of the invention to propose a solution avoiding these 

drawbacks. 

To this end, the invention relates to a video encoding method for the 
compression of an original video sequence divided into successive groups of frames (GOFs), 
said method comprising the steps of : 
20 (1) generating from the original video sequence, by means of a wavelet 

decomposition, a low resolution sequence including successive low resolution GOFs ; 
(2) performing on said low resolution sequence a low resolution decomposition, 

by means of a motion compensated spatio-temporal analysis of each 
low resolution GOF ; 

25 (3) generating from said low resolution decomposition a full resolution sequence, 

by means of an anchoring of the high frequency spatial subbands resulting from the wavelet 
decomposition to said low resolution decomposition ; 

(4) coding said fixll resolution sequence and the motion vectors generated during 

the motion compensated spatio-temporal analysis, for generating an output coded bitstream. 
30 The proposed solution is remarkable in the sense that the global structure of 

the decomposition tree in the 3DS analysis is preserved and no extra information is sent to 
correct the drift effect (only the decomposition/reconstruction mechanism is changed). If no 
motion estimation/compensation is performed at full resolution, it is a low-cost solution in 
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terms of complexity. If motion compensation is introduced in the high spatial subbands, a 
better coding efficiency is provided. 

The invention also relates to a corresponding decoding method, comprising the 

steps of : 

(1) decoding said input coded bitstream for generating a decoded full resolution 
sequence and associated decoded motion vectors ; 

(2) in said decoded full resolution sequence, separating the decoded high 
frequency spatial subbands and the decoded low resolution decomposition ; 

(3) generating from said decoded low resolution decomposition, by means of a 
motion compensated spatio-temporal synthesis, a decoded low resolution 

sequence ; 

(4) reconstructing from said decoded low resolution sequence and the decoded 
high frequency spatial subbands an output full resolution sequence corresponding to the 
original video sequence. 

The invention also relates to an encoding device and a decoding device 
provided for implementing said encoding method and said decoding method respectively. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will now be described in a more detailed manner, with reference 
to the accompanying drawings in which : 

Fig.l shows a 3D subband decomposition ; 

Fig.2 illustrates a motion-compensated temporal analysis at the lowest 

resolution ; 

Fig.3 depicts an embodiment of an encoding scheme according to the 

invention ; 

Fig.4 depicts an embodiment of a decoding scheme corresponding to the 
encoding scheme of Fig.3 ; 

Fig. 5 illustrates the reordering of the high spatial subbands (for a forward 
motion compensation) ; 

Fig.6 depicts another embodiment of an encoding scheme according to 

the invention. 



DETAILED DESCRIPTION OF THE INVENTION 
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The proposed solution (i.e. a spatial scalability with no drift in a motion 
compensated 3D subband codec) is now explained with reference to its two main steps : (a) 
motion compensation at the lowest resolution, (b) encoding the high spatial subbands. 

First in order to avoid drift at lower resolutions, Motion Compensation (MC) 
5 is applied at this level. Consequently, as illustrated in Fig.2, one first downsizes (reference d) 
the GOF using wavelet filters, and the usual 3D subband MC-decomposition scheme is then 
applied to this downsized GOF instead of the full-size GOF. In Fig.2, the temporal subbands 
(L 0 ,d, H 0 ,d) and (Li >d , Hi,d) are determined according to the well-known lifting scheme (H is 
first defined from A and B, and then L from A and H), and the dotted arrows correspond to 
10 the high-pass temporal filtering, the continuous ones to the low-pass temporal filtering, and 
the curved ones (between low frequency spatial subbands A of the frames of the sequence, 
referenced A^d, Ai >d , A2 )d , A3, dj or between low frequency temporal subbands L, referenced 
L 0 , d and Li,d) to the motion compensation (it may be noticed that a side effect of this method 
is the reduction of the amount of motion vectors to be sent in the bitstream, which saves up 
15 some bits for texture coding). Before transmitting the subbands to a tree-based entropy coder 
(for instance to a 3D-SPIHT encoder such as described for instance in the document "Low 
bit-rate scalable video coding with 3D set partitioning in hierarchical trees (3D-SPIHT)", BJ. 
Kim and al. IEEE Transactions on Circuits and Systems for Video Technology, vol. 10, n°8, 
December 2000, pp. 1374-1387), one puts the high spatial subbands that allow the 
20 reconstruction of the full resolution. The final tree structure looks very similar to that of a 3D 
subband codec such as the one described in the document "A fully scalable 3D subband 
video codec", IEEE Conference on Image Processing (ICIP2001), vol.2, pp.1017-1020, 
Thessaloniki, Greece, October 7-10, 2001, and so a tree-based entropy coder can be applied 
on it without any restriction, as described in the new encoding scheme of Fig.3, where the 
25 references are the following (for a frame of the full resolution sequence) : 
FRS : £U11 resolution sequence 
WD : wavelet decomposition 
LRS : low resolution sequence 

MC-3DSA : motion-compensated 3D subband analysis 
30 LRD : low resolution decomposition 

HS : high subbands 

U-HFSS : union of the three high frequency spatial subbands of a frame 
FR-3D-SPIHT : full resolution 3D SPIHT 
OCB : output coded bitstream. 
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The corresponding decoding scheme, depicted in Fig.4, is symmetric to this encoder (in 
Fig. 4, the additional references are the following : 

MC-3DSS : motion compensated 3D subband synthesis 

HSS : high subbands separation 
5 FRR : full resolution reconstruction). 

To enable spatial scalability, the high frequency spatial subbands just have to be cut as in the 
usual version of the 3DS codec, the decoding scheme of Fig.4 showing how to naturally 
obtain the low resolution sequence. 

Then, for coding the high spatial subbands, two main solutions are proposed, 
10 the first one without MC, and the second one with MC. 



A) Without MC 

In the first solution, the high subbands simply correspond to the high 
frequency spatial subbands of the original (full resolution) frames of the GOF in the wavelet 

15 decomposition. Those subbands allow the reconstruction at full resolution at the decoder. 
Indeed, the frames can be decoded at the low resolution. However, these frames correspond 
to the low spatial subband in the wavelet analysis of the original frames. Hence one has 
merely to put the low resolution frames and the corresponding high subbands together and 
apply a wavelet synthesis to obtain the full resolution frames. But now, where and how to put 

20 those high subbands in order to optimize the 3D-SPIHT encoder ? In a MC scheme for a 3D 
subband encoder, the low temporal subbands always look like one of the original frames of 
the GOF. As a matter of fact : 

L=-~= [A + MC(B)] . (3) 

so L looks like A. Consequently, the high spatial subband of A should be placed with the low 
25 resolution decomposition corresponding to L. This approach (reordering of the high spatial 
subband in the case of forward motion compensations) is illustrated in Fig.5, where DWT H 
denotes the high frequency wavelet filter and the coefficients cj t are multiplication 
coefficients. The way to define Cj t is described later. 

However, the motion compensation in the 3D subband structure can be either 
30 forward or backward (it has even been shown that alternate directions improve coding 
efficiency. The following algorithm, in which the notations are : 

. jt : temporal decomposition level (0 for the full frame-rate, 
jt_max for the lowest frame-rate) 
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. t : 0 for the low temporal subband, 1 for the high one 
. nf : subband index at temporal level jt 
. me_dir_descjree : a byte that describes the ME directions 
used at a given temporal level jt (the LSB describes 
the direction of the first ME/MC, 0 means "forward", 
1 means "backward"), 

makes the link between a frame GOF_index in the GOF and the spatio-temporal subband 
{jt;n;t} which resembles it most, depending on the Motion Estimation Direction Description 
Tree. 
UInt8 

STlocationToGofIndex(MEDirectionDescriptionTree me_dir_desc Jxee, UInt8 
jtjnax,UM8jt,UInt8 nf,UInt8t) 
{ 

UInt8 gof_index=0 ; 
UInt8 direction ; 
UInt8 j,n_sb ; 
UInt8 sign ; 

gofjndex = nf«jt ; 

sign= 1 ; 
nsb = nf ; 

for(j=jt-l ; j>=0 ; j~) 
{ 

direction = l«njb ; 

if(t=0) 

sign=0 ; 

direction &= me_dir_desc_tree.aui8_level[j3 ; 
direction »= n_sb ; 
if (sign) 

{ 

direction = ! direction ; 
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direction «=j ; 
gofjndex = direction; 

} 

retum(gof_index) ; 

} 

The way to define the coefficients cj t is now described (in Haar filter case). Let 
a be the coefficient used in the temporal 2-tap Haar filter. In the conventional 3D subband 
scheme, one has : 

jL = a*(A+MC 1 (B) 
|H = a*(MC(A)-B) 

If, in the present scheme, one uses c Jt = a Jt for the high spatial subbands, then 
it is still meaningful to use temporal scalability. Indeed : 

DWT L (L) = a * ( DWT L (A) + MC' 1 (DWT L (B)) ) 

DWT H (L) = c jt *( DWT H (A) ) 

= a * ( DWT H (A + UpSampl4MC A (DWT L (B))] ) 

and : 

JDWT L (H) = a * ( DWT L (B) - MC[DWT L (A)] ) 
[DWT H (H) = a*( DWT H (B) ) 

where UpSample refers to the picture upsizing using wavelet filters. For the reconstruction at 
a lower frame rate, only the low temporal subband is synthesized : 

1 DWT-'CDWTCL)] 



2*a 

| = -| * ( A + UpSamplerMC- 1 (DWT L (B))] ) 



25 



Finally, the reconstructed frames at each temporal level will tend to look like a motion- 
compensated average of the "reference" original frame and a blurred version of the other one 
(up-sampled version of the downsized frame), whereas in the current version of the 3D 
subband codec this blur is not introduced. Improving spatial scalability at the expense of 
adding blur in the temporal scalability is however a worthy step. 
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B) With MC 

As using MC in every subband does not allow a reconstruction with no drift, it 
is possible, as depicted in Fig.6, to partially use MC to construct the high spatial subbands 
(which is better in terms of coding efficiency) and still be able to reconstruct every resolution 
5 (in Fig.6, the additional references are the following : 

ME/MC : motion estimation/motion compensation 

PRE : prediction error). 
Instead of directly using the high frequency spatial subbands of the wavelet decomposition, a 
wavelet decomposition is carried out on a prediction error obtained from the MC performed 
10 on the full resolution sequence and reusing for instance the motion vectors of the low 
resolution. 

The solution is to define : 

DWT H (L) = c j t * (DWT H (A)) 
* DWT H (H) = c jt * DWT H (B - MC(A)) 

It can be noticed that the MC is only used in the high temporal subband : A is first 

1 5 reconstructed at the full resolution thanks to the low temporal subband, and then used to get 
frame B with MC thanks to H. The coefficients cj t are chosen as previously. Said MC at full 
resolution can be performed either by merely upsampling the low resolution motion vectors 
(which has the advantage of introducing no other motion vector overhead) or by refining 
these upsampled low resolution vectors (which costs some additional transmission bits but is 

20 more efficient in terms of texture coding). 

It must be understood that the present invention is not limited to the 
aforementioned embodiments, and variations and modifications may be made without 
departing from the spirit and scope of the invention. There are numerous ways of 
implementing functions of the method according to the invention by means of items of 

25 hardware or software, or both, provided that a single item of hardware or software can carries 
out several functions. It does not exclude that an assembly of items of hardware or software 
or both carry out a function, thus forming a single function without modifying the method in 
accordance with the invention. Said hardware or software items can be implemented in 
several manners, such as by means of wired electronic circuits or by means of an integrated 

30 circuit that is suitable programmed. The integrated circuit can be contained in a computer or 
in an encoder or decoder and comprise a set of instructions, contained, for example, in a 
computer programming memory or in an encoder or decoder memory and causing the 
computer or the decoder to carry out the different steps of the methods according to the 
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invention. This set of instructions may be loaded into the programming memory by reading a 
data carrier such as, for example, a disk. A service provider can also make the set of 
instructions available via a communication network such as, for example, the Internet. 
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CLAIMS: 



1 . A video encoding method for the compression of an original video sequence 
divided into successive groups of frames (GOFs), said method comprising the steps of: 

(1) generating from the original video sequence, by means of a wavelet 

decomposition, a low resolution sequence including successive low resolution GOFs ; 
5 (2) performing on said low resolution sequence a low resolution decomposition, 

by means of a motion compensated spatio-temporal analysis of each low resolution GOF ; 
(3) generating from said low resolution decomposition a full resolution sequence, 

by means of an anchoring of the high frequency spatial subbands resulting from the wavelet 
decomposition to said low resolution decomposition ; 
10 (4) coding said full resolution sequence and the motion vectors generated during 

the motion compensated spatio-temporal analysis, for generating an output coded bitstream. 

2. A method according to claim 1, in which, for each frame, said high spatial 
subbands are directly anchored to the low resolution subband that, in said spatio-temporal 

15 decomposition, looks most like said frame, depending on the motion estimation direction. 

3 . A method according to claim 1 , in which a predictive mode is used to 
construct the high spatial subbands, said high spatial subbands resulting from a second 
wavelet decomposition performed on a prediction error obtained from a motion 

20 compensation applied to the original video sequence. 

4. An encoding device for the implementation of the video encoding method 
according to anyone of claims 1 to 3. 

25 5 . A method for decoding an input bitstream coded by means of an encoding 

method according to anyone of claims 1 to 3, said decoding method comprising the 
steps of : 

(1) decoding said input coded bitstream for generating a decoded full 

resolution sequence and associated decoded motion vectors ; 
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(2) in said decoded full resolution sequence, separating the decoded high 
frequency spatial subbands and the decoded low resolution decompositipn ; 

(3) generating from said decoded low resolution decomposition, by means of 
motion compensated spatio-temporal synthesis, a decoded low resolution sequence ; 

5 (4) reconstructing from said decoded low resolution sequence and the decoded 

high frequency spatial subbands an output full resolution sequence corresponding to the 
original video sequence. 



10 



6. A decoding device for the implementation of the video decoding method 

according to claim 5. 



WO 03/063497 



PCT/IB03/00156 




WO 03/063497 PCT/IB03/00156 



2/5 




FIG.2 



FR-3D-SPlHTh >» 



LRD 
LHHFSS 



MC-3DSS 




f 



HSS 



FRR 



FRS 



RG.A 



WO 03/063497 



PCT/IB03/00156 



3/5 




CD 



WO 03/063497 



PCT/TO03/00156 




WO 03/063497 



PCTYIB03/00156 




CO 



INTERNATIONAL SEARCH REPORT 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 H04N7/26 



According to International Paient Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 H04N 



Documentation searched other than minimum documentation to the extent that such documents are included In the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

EPO-Internal 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category • Citation of document, with Indication, where appropriate, of the relevant passages 



Relevant to claim No. 



BOTTREAU V ET AL: "A FULLY SCALABLE 3D 
SUB BAND VIDEO CODEC" 

PROCEEDINGS 2001 INTERNATIONAL CONFERENCE 

ON IMAGE PROCESSING. ICIP 2001. 

THESSALONIKI , GREECE, OCT. 7 - 10, 2001, 

INTERNATIONAL CONFERENCE ON IMAGE 

PROCESSING, NEW YORK, NY: IEEE, US, 

vol . 2 OF 3. C0NF. 8, 

7 October 2001 (2001-10-07), pages 

1017-1020, XP001045747 

ISBN: 0-7803-6725-1 

cited 1n the application 

the whole document 

-/-- 



1-6 



m 



Further documents are listed In the continuation of box C. 



Patent family members are listed in annex. 



• Special categories of cited documents : 

'A' document defining the general state of the art which Is not 
considered to be of particular relevance 

'E' earlier document but published on or after the International 
filing date 

■L' document which may throw doubts on priority claim(s) or 
which Is cited to establish the publication date of another 
citation or other special reason (as specified) 

'O' document referring to an oral disclosure, use, exhibition or 
other means 

■P" document published prior to the international filing date but 
later than the priority date claimed 



"T later document published after the international filing date 
or priority date and not In conflict with the application but 
cited to understand the principle or theory underlying the 
Invention 

■X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
Involve an inventive step when the document is taken atone 

*Y* document of particular relevance; the claimed invention 

cannot be considered to involve an Inventive step when the 
document Is combined with one or more other such docu- 
ments, such combination being obvious to a person stdBed 
In the art 

document member of the same patent family 



Date of the actual completion of the International search 



23 May 2003 



Date of mailing of the international search report 



02/06/2003 



Name and mailing address of the ISA 

European Patent Office, P.B. 5816 Patentlaan 2 
NL-2280 HV Rijswijk 
Tel. (431-70) 340-2040, TX. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Authorized officer 



Schoeyer, M 



Form FCT/ISA/210 (second sheet) (July 1992) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



PCT/IB 03/00156 



C(Contlnuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 


Category • 


Citation or document, with indicatlon.where appropriate, of the relevant passages 


Relevant to claim No. 


X 


WU F ET AL: "DCT-predict1on based 
progressive fine granularity scalable 
coding" 

PROCEEDINGS, INTERNATIONAL CONFERENCE ON 
IMAGE PROCESSING, XX, XX, 
10 September 2000 (2000-09-10), pages 
556-559, XP002165186 
the whole document 


1-6 


A 


US 2002/006164 Al (FELTS BORIS) 
17 January 2002 (2002-01-17) 
abstract 

page 1, paragraph 1 -page 3, paragraph 46 


1-6 



Form PCT/ISA/210 (continuation of second ahoet) (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

patent family members 



PCT/IB 03/00156 



Patent document 
cited In search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



US 2002006164 Al 



17-01-2002 CN 
WO 
EP 



1381142 T 
0189226 Al 
1287701 Al 



20-11-2002 
22-11-2001 
05-03-2003 



Form PCT/ISA/210 (patent family annex) (July 1998) 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 



Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 



□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 



IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




COLOR OR BLACK AND WHITE PHOTOGRAPHS 



